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ABSTRACT 



We contrast the two approaches to "classical" signature change used by Hayward with the 
one used by us (Hellaby and Dray). There is (as yet) no rigorous derivation of appropriate 
distributional field equations. Hayward's distributional approach is based on a postulated 
modified form of the field equations. We make an alternative postulate. We point out an 
important difference between two possible philosophies of signature change — ours is strictly 
classical, while Hayward's Lagrangian approach adopts what amounts to an imaginary proper 
"time" on one side of the signature change, as is explicitly done in quantum cosmology. We 
also explain why we chose to use the Darmois-Israel type junction conditions, rather than the 
Lichnerowicz type junction conditions favoured by Hayward. We show that the difference in 
results is entirely explained by the difference in philosophy (imaginary versus real Euclidean 
"time"), and not by the difference in approach to junction conditions (Lichnerowicz with specific 
coordinates versus Darmois with general coordinates). 
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DIFFERENCES IN APPROACH 



Meaning of 'classical' 

The most fundamental difference between the approach of Hayward and that of Hellaby 
and Dray (us) is in the treatment of the region of Euclidean signature. We (Hellaby and Dray) 
[1,2] adopt a strictly classical view of this region, as do Ellis et al [3-7] and Dray, Manogue and 
Tucker [8-11], whereas Hayward [12-14], in his Lagrangian approach [12], works with squared 
momenta that change sign, which is equivalent to retaining the imaginary proper "time" and 
momentum fields that appear in the quantum cosmological treatment [see e.g. 17,18], obtaining 
a classical version of the quantum cosmological result. In an interesting comparison, Kerner 
and Martin [20-21] have done a quantum cosmological version of the strictly classical signature 
change, retaining the real "time" of the classical treatment of Euclidean signature, so that the 
wave function is oscillatory in the Euclidean region as well as the Lorentzian region. It is not 
surprising that these four approaches (classical or quantum cosmology with real or imaginary 
Euclidean "time" ) all lead to somewhat different results. None of them can claim to have any 
experimental or observational support. 

Appropriate form of the EFEs 

The Einstein Field equations were derived assuming a constant signature metric. For sig- 
nature change, the most basic question to be addressed is that of the applicable form of the 
Einstein equations at a change of signature. Various postulates have been made and will be 
discussed below. The differences between them can be reduced to the question of whether the 
Euclidean proper "time" is real or imaginary. The applicable form of the conservation equations 
is similarly open to question. 

Distributional EFEs 

Formal calculation of the Einstein tensor from a step function metric results in products 
of distributions, so another approach is required to obtain a well defined distributional form 
of the Einstein field equations (EFEs). Hayward has shown [12,14] that the usual form of the 
EFEs for a lapse function metric may be interpreted distributionally, but the derivation of these 
equations requires certain (unstated) smoothness assumptions. This does not in itself support 
those assumptions or that form of the distributional EFEs. We postulate an alternative form 
for the distributional EFEs based on more physical (and less restrictive) assumptions. 

The approach of Kossowski and Kriele [15-16] also makes restrictive smoothness assump- 
tions. Their comparison of two "classical" approaches investigates under what conditions they 
are equivalent, showing that in general they are not; and in each case under what conditions 
the Einstein tensor is bounded and/or continuous. However, they require the full 4-dimensional 
matter tensor to be not only bounded but also smooth through S. This means their proposal 
(a) is not what is done by us and by Ellis et al [3], so that their remark 2 criticising [3] and 
their very restrictive results from theorems 9 and 10 are not applicable. 

Choice of junction conditions 

The standard view for realistic Lorentzian solutions of Einstein's field equations is that 
the metric should be C 1 (on one or more 3-surfaces) and piecewise C 3 (everywhere else). The 
surface of the moon is an example of a timelike surface dividing regions of different density. The 
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time of decoupling in the early universe is an example of a spacelike surface dividing regions of 
different pressure. In other words, some components of the Einstein/matter tensor have jumps 
at junctions. How should a surface of signature change be treated? It is entirely reasonable 
to permit similar or geater discontinuities at a signature change — a change which is far more 
fundamental than a change in the equation of state. 

In the Lagrangian approach [12], the arguments for the vanishing of the momenta are similar 
to those used in quantum cosmology [17,18] for the vanishing of the extrinsic curvature. To 
satisfy the Einstein field equations (EFEs) and conservation equations, through the signature 
change surface E, implicitly requires the metric and the matter fields (in this case a scalar field 
0) to be C 3 . In contrast, the Lichnerowicz junction conditions for the case of no signature 
change (which are equivalent to the constant signature Darmois junction conditions) say that 
if you can find admissible coordinates spanning E in which the 4- metric is C 1 , then you have 
a good matching. At a signature change, this continuity (whether C 1 or C 3 ) also requires 
the choice of a particular class of coordinates — lapse-shift coordinates with a squared lapse 
function, such as N(t) = t. We point out below that this unusually high degree of continuity 
can only be achieved in a limited sense, owing to the fact that the metric is degenerate at E, 
meaning the inverse metric is divergent. 

The alternative choice of the Darmois junction conditions is attractive because they are 
expressed in terms of invariants of the signature change surface E — the intrinsic (3-d) metric 
jjij and the extrinsic curvature K^. In the familiar constant signature case, if the one-sided 
limits of these two on one side of E (i.e. in manifold with boundary M + ) equal the one- 
sided limits on the other side (M~), then there is a good matching, and it is possible to find 
coordinates spanning E which make the metric and the inverse metric at least C 1 . The (one- 
sided limits of these) invariants do not depend on the choice of coordinates in the enveloping 
manifolds M~ and M + . In particular, for a spacelike surface, they do not depend on the 
choice of time or lapse function. For the case of a change of signature at E, careful analysis [1] 
showed that no modification of these conditions, or any special coordinate choice is required. 
The Darmois conditions are blind to the change in signature, just as the EFEs are blind to 
the metric signature. Of course it is no longer possible to find coordinates that make both the 
metric and its inverse C 1 , but, as above, the covariant metric can be made arbitrarily smooth. 
Thus all the main results in [1] are coordinate invariant. Supplementary junction conditions 
on the equation of state or the matter fields may also be imposed, and such conditions were 
given in [1] for an electromagnetic field, as well as the consequences of the EFEs at a signature 
change in an arbitrary electrovac model, in equations (78)-(80). 

Standard conservation can be shown to hold as a consequence of the constant signature Dar- 
mois (or Lichnerowicz) conditions. At a signature change, the same conditions lead to extra 
jumps in physical quantities, such as the density, and a surface term in the derived conservation 
law. Conservation is not abandoned, but it has to be modified. Stronger conditions than those 
in [12] are required to restore the conservation law to its familiar form, for the strictly classical 
case. 

THE CASE OF A HOMOGENEOUS SCALAR FIELD 

The Field Equations for a squared-lapse metric 

Consider the case of a homogenous isotropic universe with scale factor a(t) and squared lapse 
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function N(t), filled with a scalar field <fi(t) with arbitrary potential V(4>). The line element 




ds 2 = -N(t)dt 2 + hdy'dy 3 = -N(t)dt 2 + a 2 (t) ( 7 dr \ - + r 2 (d9 2 + sin 2 6d X 2 )) (1) 

[ (1 — At 2 ) J 

has A" C° and N(t = 0) = on the signature change surface E. The and rr EFEs Gp = nTp, 
evaluate to [24] 

where ' = d/dt, the 99 and XX equations adding nothing new. The covariant and contravariant 
forms re-scale the tt but not the rr equations by factors of N. The Klein-Gordon equation is 
implied by the EFEs. 

Strictly speaking, these equations are not well defined on E except in a limiting sense as 
N — > 0, and for k and V non-zero, G a = kT 11 is divergent in this coordinate system, regardless 
of the behaviour of a 1 and near iV = 0. Similarly, the conservation equation V Q G^ or V Q G a ^ 
is badly defined at iV = 0. The latter form is especially problematic, since G tt is divergent for 
k 0. One could therefore question whether the EFEs and the conservation equation are well 
defined on E. A reasonable possibility is to evaluate these equations away from = and 
rescale by factors of N, to obtain equations which are well defined at A^ = 0. The 'second order 
equations' of [12] and equations (2) of [14] are essentially re-scalings by factors of N, of the 
Einstein, Raychaudhury, Klein-Gordon, and/or Friedmann equations. Note that the re-scaling 
itself is only well-defined in the limit as t — > 0, not at t = 0, itself. 

The Darmois approach with real "time" 

In this approach, we postulate that two metrics of different signature, each of which solves 
the respective constant signature form of the EFEs, may be pieced together so that the Darmois 
junction conditions are satisfied at the signature change surface E. Even if the EFEs are well 
defined at E in each solution, the resulting matter fields are not necessarily smooth when put 
together. This is true with or without signature change. 

We patch together the Lorentzian and Euclidean signature forms of the above metric (]T|) 
(M~ and M + ), identifying the two t — surfaces in each to form the junction surface E. Using 
the notation of [1], the strictly classical treatment requires a real proper "time" 

+1, t < 0, A^<0 Euclidean signature metric 
dr = V-eN dt , e = { (4) 

—1, t > 0, A^>0 Lorentzian signature metric 

The choice of the Darmois approach, in which all quantities are well defined without re-scaling, 
in the limiting sense given in [1], is made largely because of this extra measure of discontinuity, 
and partly because of the above-mentioned singularities in the metric and Einstein tensors. In 
M ± we calculate 

% = diag(> rr , gee, g xx ) = a 2 diag j^y , r 2 , r 2 sin 2 d\ (5) 
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,0,0,0 



K i3 



K 



{V-eN ,0,0,0) , lj a = e 



= di&g(K rr , K ee , K xx ) = 



aa 



-diag 



2 2 • 'A n 

r , r sin 9 



kr 2 ) 



3a' 



a^/-eN 
6k 



= 3 



' d(lna) 








T^l^lt, = eT ° = — 



U\2 



N 



+ V 



1 

2 



(1 - Ax 2 ) 



T 01 = T* = Tj?l lt e'[ = 



(6) 
(7) 
(8) 

(9) 

(10) 

(11) 
(12) 

(13) 
(14) 



where n a and l a are chosen to point from M~ into M + . The tilde indicates a quantity evaluated 
in geodesic normal coordinates based on E, used for the purpose of easy interpretation, and for 
compact notation indicating invariants. With the notation 



IE ' 



Lim 



(15) 



for the jump in the value of Z across E, and choosing to keep r, 9 and <fi continuous through 
E, the first Darmois conditions = give the obvious relations 



k+ - k~ = [k] = , 
and the second Darmois conditions [K^] = give 



a i - a s = N = ° 



(16) 



(17) 



where r is the comoving proper "time" . Appropriate supplementary junction conditions for the 
scalar field are 



a' 


+ 




a' 


da 
dr 


+ 


da 




da 


V-eN 


e V-eN 


s 


V-eN _ 


E 


dr' 


E 


dr 



[0] = , [n a d a <P] 



d(f> 




\ ' 1 


dr 







= , 



[V] = o , 



dV 



(18) 



From eqs (47)-(50) of [1], the various surface effects in the conservation law are then 



E(l a ) = [G 00 ] = -ft 
E(el) = [G 01 ] = 2( 3 V J X lj - V,K) 
E(n a ) = [G° ] = K 2 - KijlC* 



-% = k[T 00 ] = -kVz 

a E 

= k[T 01 ] = 



-eN 



= k[T°] = 



(19) 
(20) 
(21) 

(22) 
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Assuming that £ is not already a singular surface in one or both of M~ or M + separately, 
these quantities are all well defined in the limit as iV —>■ without re-scaling. From the EFEs 
we obtain the restrictions 

6 \dr)^ ~ { dr J s ' 2 a| 1 j 

and while (da/dr)^ = and (d<j)/dr)-£ = are possibilites, they are not demanded. In any case 
it is not possible to remove the surface effect E{l a ) = — 3 i?| s = — (6A;/a|) for A; 7^ 0, in this 
scenario. 

In summary, we are given a and of the form 

a = a + (t)B(t) + a_(t)B(-t) , (f> = <f>+(t)B(t) + <f>-(t)B(-t) (24) 

such that a±, <j>± G C 3 separately solve the constant signature EFEs, as well as satisfying the 
Darmois and Klein-Gordon field junction conditions (|T6"|)-(|l8"|) and signature change condition 
fl23|) . We also assume 

a'+ <j)'+ , , /2a'| a'+iV'N „ . , 

lim === =■ , lim t^= - and lim — all exist. (25) 

t-o v^iiV t-o v^iiV *-o ^ AT AT 2 / v y 

It follows from (0) and (0) that [a] = = [0], so 

a' = a' + (t)e(t) + a'_(t)e(-t) + [a]6(t) = a' + (t)Q(t) + a'_(t)<d(-t) (26) 

0' = <p' + (t)&(t) + #_(*)e(-f) + M*(t) = 0' + (t)e(t) + #_(*)©(-*) (27) 

a" = al(t)e(t) + a'L(t)e(-t) + [a')5(t) (28) 



While condition (17), = [d] = [a']/y\N\ , certainly ensures [a'] = 0, we must verify that it is 
sufficient to eliminate distributions when we insert this into (0) and (|3]). For example, putting 
' = d/dr, 

[a]=0 = [o] d = a+Q(t) + 2_6(-t) 

/ iV'al a" \ _ , (N'a' a" 



V 2iV 2 A r / 
/ N'a' a"\ 




2N 2 N 



(29) 



and there is thus no distribution in a"/N. Conversely, if ((N'a')/(2N 2 ) — a"/N) is free of dis- 
tributions, then ed is. Therefore all terms in the EFEs are bounded, and we have a solution of 
(fj) and (Q), as well as the Klein-Gordon equation. We stress that the Darmois conditions are 
imposed on a(t), 0(£) and V(4>) before calculating and from them. Our assumptions are 
borne out by solutions to (@)-(^D which will be exhibited elsewhere [27]. A rigorous treatment 
of distributions in the presence of signature change will appear in [26] . 

The distributional EFEs 

In the distributional approach, we must postulate a form of the EFEs that makes sense 
when N(t) becomes the the step function — G(— t) + 0(i) = e(t) and N' becomes 25, where 
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0(t) is the Heaviside function and S(t) is the Dirac delta function. In particular products of 
distributions must be avoided. Then t = r and e = —e. There is currently no known way to 
derive such equations from first principles. 

Postulate A. The re-scaled EFEs, A/"x(0) and iV 2 x(|3|), are such a form if it is assumed that 
a, a', (ft, 4>', V(<f>) (and dV/d<p) are at least C°, and a" (and </>") may contain a step function, 
but no distribution. It then follows [12,14] that there is a distributional term in iV 2 x G r r , but 
not in iV 2 x /tT r r , unless a' s = 0. Given these assumptions with this postulated form of the 
distributional EFEs, a' s = is required. 

However, inspection of a" in (|31~D and in solutions to (@)-(§) [27], with N — > e and N' — > 25, 
suggests that this assumption may be too strong, and makes it questionable that a" — > (d 2 a/dr 2 ) 
without any extra distributional terms. In other words, the distributional term that appears 
in 7V 2 x(3) when N' — > 5 is only present because (|2])-(|3]) were derived and re-scaled using the 
differentiable structure of a continuous N, which is not appropriate for a step function e. 

Postulate B. Recognising that this substitution effectively replaces lapse "time" with proper 
"time" ,t—*T, which affects the functional form of a, we consider it safer to transform from t to 
t 7^ t, rather than substituting functions with distributions. We write ' = d/dr," = d 2 jdr 2 and 
we continue to assume the Darmois and scalar field junction conditions, (p76|) - (p~8|) , to eliminate 
any distributions in the following transformations: 



T+ = fiy/N{t)dt, t>0 
r(t) = { (30) 



T- = -J?y/-N(t)dt, t<0 

f VNd, t>0 \ Nii+^Ld, t>0 

a = a(r(t)) , a' = I , a" = I (31) 

{ a , t < [ -Nd - , t < 

r s/n $ , t > o 

= 0(r(t)) , 0' = (32) 

{ ^N<j), t < 

Using these expressions in d26|)-(p8|) and substituting the result into (0) and (|3|) leads to: 

G < = "I i k - e <°> 2 } = -^{-fW 2 + V} = < (33) 
r £ = - H ^_ v} = (tJ? (34) 



where we have again used (p9|). We note that this form of the equations does not contain iV or 
N', nor does it need rescaling, so it seems more natural for a proper time solution. We have 
assumed that a, a, <p, <p, V ((/>), (and dV/dip) are at least C°, so that a (and 4>) may contain a 
step function but no distribution, and we postulate (|33|) and ( p4[) for the field equations when 
N — > e. All of the solutions in [27] satisfy these EFEs. In particular, = is not required. 

The Lagrangian approach with imaginary "time" 

We postulate a Lagrangian approach in which the squared momenta are required to be 
smooth. The action for this case, in terms of the gravitational and Klein-Gordon Lagrangian 
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densities can be written as 

S = J(C G + 2KC KG )d 4 x (35) 



Cg = J=~gR = JW{ i R + K ij KV-K 2 ) (36) 



Ckg = ~^y~ {(9^)9^(9 A) + V(<I>)} (37) 

where R and g are the Ricci scalar and the metric determinant, and ?/ = a 6 r 4 sin 2 9/(1 — kr 2 ) 
is the determinant of the 3- metric ^ of (|l|). Leaving the 4- volume element y/—g d 4 x in this 
form leads to an effective imaginary proper "time" 



da = VN dt (38) 

in the Euclidean region, N < 0. This leads to an imaginary unit "normal" n a = dx a /da = 
(1/y/N, 0, 0, 0) for N < which still obeys n a n a = —1. In this way this region is not Euclidean 
in the classical sense. It also means the extrinsic curvature and the momenta change from 
imaginary to real across iV = 0, as noted by Hayward [12], 

* = § = n"^= ****** (39) 
d<P y/N{l - kr 2 ) 

K tj = n a V a hj = ^f ( 40 ) 



. N 
BC r- 

P) = hkgjjT = \H)( K j- K V) , K = K Z (41) 



Hk 

2a 2 a'r 2 sin9 



ay/N(l - kr 2 ) 3 
and the squared momenta change sign, 



(42) 



2 M'Wr* sin 2 9 , s 

(P% = p)p{ = K)K{-2KKl + K 2 5l (44) 
4(o') 2 aV sin 2 tf , 



N(l - kr 2 ) 



SI (45) 



implying if)*, (p*)^ , ip, p u and must all be zero at the change surface assuming the momenta 
are required to be smooth and a, <f), <$' and a' all real. 

The Lagrangian approach with real "time" 

One could argue that it is more reasonable to define the action using a real volume element 
\g\ d A x 



S = J(C G + 2KC KG )d 4 x (46) 



C G = J\g\R=Jg\N\(*R + K ij KV-K 2 ) (47) 



Ckg = -^r- {(<9^W) + V(<j))} (48) 



and this form leads to a real proper "time" everywhere 

dr = J\N\ dt (49) 



so that n a = dx a /dr = (l/y\N\ ,0,0,0) is real and genuinely spacelike n a n a = +1 in the 
Euclidean region. The momenta now stay real 

<j)'a 3 r 2 sin 9 , , 

i) = =■ (50) 

y/\N\{l -kr 2 ) 

K i3 = (51) 
N 



P) = ~ 2aVr2sing 5) (52) 



ay/\N\(l - kr 

and the squared momenta don't change sign, 

_ (0QWsin 2 fl 
^ " \N\(l-kr*) {b6) 
. 4(a') 2 a 4 r 4 sm 2 6 ■ 

Thus the above requirements of smoothness and reality lead to the conditions that (a') 2 /\N\ 
and (<j)') 2 /\N\ be smooth, in agreement with the real "time" Darmois based approach. 

The Darmois approach with imaginary "time" 

We now apply the Darmois based approach to the quantum cosmology inspired model of 
signature change with an imaginary proper "time" in the Euclidean region, using the same 
metric. We list quantities that are different from above. 

da = \fN dt = 4^1 dr (55) 

n a = 0,0, (A , h a h a = -l (56) 



.00 



1« = (VN, 0,0,0) , ("i° = -l (57) 

*« - w Aias {j^) yyAn2g ) (58) 

where the double tilde means that the geodesic normal proper "time" coordinate a is imaginary 
where N < 0. Applying the Darmois junction conditions as before, and again requiring a, 0, 
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a' and 0' to be real, now gives 



[k] = , [a] = , 



(ia 
d<r 



(62) 



which leads to 



k + = AT 



da 
''aV 



da 
aV 







(63) 



by the usual argument [18]. The supplementary scalar field junction conditions 

~dV~ 



[0] = , [n a <9 Q 0] 
lead to 



da 







[V] = , 





+ 






d0 


+ 


d0 


d0 




d0 


s 




£ 


d7 



d0 



= = 







0' 



(64) 



(65) 



Clearly <p' and a' must go to zero faster than V iV . These requirements, derived using the imag- 
inary Darmois approach, are in agreement with those of the imaginary Lagrangian approach 
and with quantum cosmo logical type signature change. This approach guarantees that G^ v and 
are smooth and finite, and conservation is observed, since the surface effects are all zero 

too 



E(l a ) = — 
a E 





1 


m 2 ] 


. N 


~ 2 


N 



, E(ei) = , E(h a ) = -E(l a ) = , f?(e?) = 

(66) 



Of course G tt and T a are still divergent. 



The EFEs with real "time" 

When the EFEs (2)-(3) are calculated from (1), it is not necessary to specify the real or 
imaginary nature of the "time" direction, nor is it evident upon inspection. This is because a 
smooth signature changing metric does not contain information on whether tangent vectors to 
geodesies in the "time" direction change character across S. 

In order to implement strictly classical signature change, we must include the fact that n a 
changes from a (real) space-like direction to a (real) time-like direction; n^g^ v n u = —e. Taking 
the derivative of this along n a gives [5] 

n»n v n x V x9ilv + 2n^X = 26(t) , a v = n x V x n u (67) 

So on S, either a u is not orthogonal to n u , or V\g^ u is not zero, and the latter seems more 
likely given that signature change is a metric phenomenon. In any case, this introduces an 
extra degree of discontinuity not revealed by calculating the Einstein equations directly from 
the metric, so that we have to write down the results on either side of E, and then carefully 
examine the junction. The fundamental role played by the unit normal was stressed in [2,6]. 
Also, because of the metric degeneracy, a classical relativist suspects the t coordinate may be 
badly behaved, just as the Schwarzschild time T is, in the exterior Schwarzschild metric at 
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R = 2M. Thus, with the EFEs (@)-@ calculated directly from the metric (1) we project into 
n a 's frame, obtaining 

& = ^ = ^l a n^^-k-^ = |j-M)!-y} (68) 

e 2 G 00 = G 00 = (?*l a l p = KT^Up-* I {(e)* + = ||Mp! + ( £ )y|(69) 

Assuming that the values of k, a, a', N, N', (ft', and V are continuous through E suggests the 
signature change conditions ( p3|) in order that the Einstein equations hold on both sides of E. 
However, going from Lorentzian (e = +1, AT > 0) to Euclidean (e = —1, AT < 0) regions, G 00 , 
Too, G 00 , and T 00 clearly jump by — (6fc/a|) = — kVe — as in (|T9| ) — and this jump cannot 
be removed for k ^ 0. Since the quantities in ( p8[) are all finite (unless E is already singular 
in M^), and since the sign of A" changes across E, it follows that G®, and Tq also jump by 
6(d{\aa)/dr)l = «(d0/dr)| — as in (H). 

Similarly, projecting the momenta fl39|) and fl42| ) into the "time"-like frame defined by n a , 
gives us (|50D and (^) again. 



CONCLUSION 

We have just demonstrated that (a) the results of the strictly classical Darmois approach 
to signature change remain the same when squared-lapse coordinates are used; (b) when the 
Darmois approach is applied to the imaginary "time" model of signature change, it leads to 
Hayward's results; (c) when the Lagrangian approach is applied to the strictly classical model 
of signature change, our results are recovered. 

There is a basic conceptual problem with this imaginary "time" . Coordinate invariance and 
local Lorentz invariance are highly desireable features of classical Relativity. Can they be shown 
to hold in general for a metric structure with positive definite signature that admits vectors 
with negative magnitudes, in some physically meaningful sense? Are the Einstein equations well 
motivated? How do we know what is 'physically meaningful' in spaces with these properties? 
These issues need scrutiny. An interesting possible approach has been given in [22-23], in which 
the "time" direction is complex, and the Wick angle is treated as a dynamical degree of freedom. 

Clearly the difference in results is entirely due to a different philosophy of signature change, 
and not at all due to the different approaches to junction conditions, or to the choice of coor- 
dinates. Strictly classical signature change does produce surface terms in the conservation law, 
but the imaginary time signature change has a Euclidean region with time-like features. 

Hayward's approach focusses on maximising the smoothness of the covariant metric, the 
covariant Einstein tensor and the matter fields through E in a particular type of coordinate 
system. Our approach focusses on the smoothness of the geometry and what kind of jumps in 
physical measurables, i.e. scalar invariants, are allowed across E. In his Lagrangian approach, 
Hayward effectively chooses imaginary "time", and his method then obtains well-defined re- 
scaled EFEs through the signature change, at the cost of not-well-defined G^ u and T^ u and 
divergent g^, G^ v and T^ v '. We (Hellaby and Dray) choose real "time" and our method obtains 
well-defined, finite, coordinate-invariant expressions for the jumps in certain physical quantities 
at the signature change surface, at the cost of finite but discontinuous g^ u and g^ u . Even in 
the absence of surface terms, bounded jumps in G^ u and T^ u are standard at boundaries. 
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We leave the reader to weigh the two philosophies of signature change, and the problems 
associated with each. 

The EFEs for a step function metric are badly defined. There is (as yet) no unique pre- 
scription for obtaining well defined distributional field equations. Hayward's postulated form 
does lead to his results, but requires a re-scaling of the EFEs. We postulate a form that doesn't 
require re-scaling and leads to our results. In the absence of a derivation of the distributional, 
signature-changing EFEs from first principles, one should be careful not to claim that a par- 
ticular form of these equations is "the" EFEs. Rather, one must investigate and compare the 
properties of alternative definitions. We emphasize that in addition to being unrescaled, (|33| ) 
and (|34l) are equivalent to (0) and (f|) for continuous N, whereas the distributional form of 
N 2 x (HP contains an extra distributional term that is not removed by the associated coordinate 
tranformation. 

We regard Hayward's approaches to signature change as reasonable and interesting, but not 
the only possibilities. 
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